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Abstract 

We study radiative px-broadening of high energy quarks passing through hot and cold QCD 
^^ ' matter. With L the length of the matter and Iq the size of constituents of the matter we find (pj^) 

^N . has both double logarithmic terms, In (L/Iq), and single logarithmic terms, ln(L//o), coming from 

gluon radiation induced by the matter. We use a (slight) extension of a formalism developed by 
^ I B. Zakharov for studying energy loss, a formalism which, for much of our calculation, reduces to a 

simple dipole scattering analysis. We estimate the radiative contribution to be a sizable correction 
^SJ ' to the nonradiative value of {p^)- We also carry out a resummation of the double logarithmic 

terms that we find, and we briefly discuss running coupling effects which appear here in a rather 

r^ . unusual way. 
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I. INTRODUCTION 

High energy quarks and gluons passing through hot and cold nuclear matter have their 
transverse momentum distributions broadened due to multiple, inelastic, scattering with 
the constituents of the nuclear matter [l|, |2| . This transverse momentum broadening can be 
nicely expressed in terms of the elastic scattering of a quark or gluon dipole in the nuclear 
matter as given in Eq. ([1]) below. 

In fact the relationship between, say, quark transverse momentum broadening and the 
elastic scattering of a quark-antiquark dipole given in Eq. ([T]) appears to be a rather gen- 
eral relationship in gauge theories [31,^]. In this paper we will exploit this relationship to 
evaluate p^-broadening where the broadening may come from gluon radiation as well as 
from multiple scattering. More precisely, p_L-broadening from inelastic multiple scatterings 
is given as p\ = qL where L is the length of the material and q is the transport coefficient 
corresponding to the hard parton, quark or gluon, passing through the medium, hot or cold 
QCD matter. In this paper we evaluate next to leading order corrections which are enhanced 
by single and double logarithms, that is radiative corrections of size {asNcqL/'K)\ri'{L/lQ) 
and {asN^qL/ir) ln(L//o) where Iq is the "size" of the scattereres in the medium. We shall 
see explicitly that our answer is independent of the exact definition of Iq. 

To a large extent this paper is motivated by an earlier paper by one of us (BW) who 
dealt with the same problem in Ref. |5|. That paper formulated the next to leading order 
Px-broadening problem in terms of the BDMPS formalism |2|, |6|, |7| developed for calculating 
radiative energy loss of high energy partons passing through nuclear matter. However, the 
BDMPS formalism has a shortcoming, namely the single scattering term with the medium 
is absent. For radiative energy loss this is not a difficulty, at least for a sufficiently large 
medium, since the dominant contribution comes from multiple scattering. But in the present 
context the single scattering term is crucial. There is a second formalism, developed by 
B. G. Zakharov [8l,l9l] for calculating radiative energy loss. While the formalism of Zakharov, 
in contrast to BDMPS, does not explicitly eliminate the single scattering term the final 
answers for radiative energy loss are identical in the two formalisms. However, for radiative 
Px-broadening this is not the case and as a result Zakharov formalism seems better suited 
to the broadening problem. We could get our result Eq. (H2]) by taking the calculation of 
Ref. |5| and adding the single scattering term missing in Ref. [5| . We have done this but we 
believe our present discussion to be simpler. 

Zakharov formalism is a dipole-like formalism with the difference being that the quark- 
ant iquark-gluon system, dominant in radiative energy loss or radiative p^-broadening, is not 
frozen in impact parameter over the times of the scatterings as in the case in normal high 
energy scattering uses of the dipole formalism Uj, |l2|]. However, in the single scattering 



approximation in our current application a Zakharov type formalism is the same as the 
normal dipole formalism. As we shall see below all the double logarithmic terms and some 
of the single logarithmic terms come from single scattering and must be given by a simple 
dipole evaluation. 

It is not hard to see why the single scattering plays such a prominent role in the present 
calculation. In the dipole picture nonradiative p^-broadening is given as the Fourier trans- 
form of the elastic scattering of a dipole of transverse size x± with the medium as given 
in Eq. ([1]). The dipole size which dominates the integral in Eq. ([1]) is x_l ~ 1/Qs where 
Ql = qL so that in terms of the dipole formalism only a few elastic scatterings, over the 
length L of the medium, are used to determine the typical transverse momentum given by 



the many inelastic scatterings of the quark passing through the medium. If we now add a 
radiative gluon to the dipole the gluon can be characterized by its hfetime and the transverse 
size {B±) that it makes with the original quark-antiquark dipole of size x±. In order to get 
a double logarithmic contribution it is necessary to have logarithmic integrals over both the 
lifetime of the gluon fluctuation and -B^. It is clear that Bj_ must be much larger than x±, 
because the quark-antiquark dipole only has, typically, a single scattering over a time L and 
now we wish to have the quark-antiquark-gluon system to interact with the medium over 
a time, or distance, short compared to L. However it is known that when B± 3> x± the 
probability of a quark-antiquark dipole to emit a gluon goes as x\dBj_/Bj_. This is not 
a logarithmic integral, however, it becomes logarithmic if and only if one includes a single 
scattering, which gives an additional factor of Bj^. 

We are now in a position to indicate which kinematic regions give single and double 
logarithmic contributions in terms of the regions shown in Fig. [2] where, for technical conve- 
nience, we have chosen to characterize the radiative gluon by its lifetime, t, and its energy, 
u. As indicated in the figure the lifetime can be between zero and L while the energy of 
the gluon is between and gL^. The interior of the region between lines (a), [b) and (c) 
gives the double logarithmic contribution to the radiative contribution to {p\) of the quark. 
Single logarithmic contribution comes from regions close to the lines (a), (b) and (c). The 
line (6) is the transition between single scattering dominance, below (6), and multiple elastic 
scattering where the S in Eq. ([1]) becomes small, above (6). The region close to line (a) is 
a single scattering region where the transverse distance between the gluon and either the 
quark or antiquark is comparable to the quark-antiquark separation x±. The region around 
line (c) is the region where the lifetime of the gluon is on the order of the proton's size, for 
cold matter. In this region the transport coefficient can no longer be treated as a constant 
since the gluon distribution included in q is evaluated at large x-values. This region is, 
however, dominated by single scattering with the medium. 

We should comment on some of the limitations and subtleties of our calculation, (a) 
We have done our calculation in a context of taking the transport coefficient q to be scale 
independent. This is not necessarily the best choice, so in Sec. IVIII Bl we do a running 
coupling calculation of the double logarithmic radiative terms. Curiously the double loga- 
rithmic contribution to {p\) turns out to be exactly a factor of two greater in the running 
coupling case as compared to fixed coupling and, remarkably, no In In terms appear. We do 
not understand how to do single logarithmic terms in the present running coupling context. 
(b) In all the regions in which there are important contributions it is only near the boundary 
(a) in Fig. [2] that the emitted gluon transerverse coordinate gets as close to either the quark 
or antiquark as the size of the quark-antiquark separation. Thus except for the boundary 
(a) calculation given in Sec. IVIBI we can define a mean p^, {p\), using Eq. (TT71) . However 
near the boundary (a) there is an x\ \nx\ dependence for S{x±) so that the mean oi p\ no 
longer makes sense. This corresponds to a p^ behavior for large p\ in dN/(P'p^_. An almost 
identical issue has been discussed earlier in Ref. |2[ and Ref. [l3| where it was suggested to 
replace the idea of a mean p\ by that of a typical p^. This amounts to taking the V^ in 
Eq. ( IT71) on all terms other than the In x\ and at the end setting x\ to its normal scale, in 
our case x\ ~ 4/Q^. In the calculation given in Sec. IVIBI we follow this prescription and 
leave the answer, in Eq. ( pTl) and Eq. (H2l) . in terms of x\. In our final estimate we take 
x]_ = 4:/Ql to get Eq. fj43l) . (c) The double logarithmic terms can be resummed, and this 
is done in Sec. IVIII Al However, it seems that one needs media lengths rather large for this 
resummation to be useful. For lengths of about 5fm the first leading double logarithmic 



term gives the dominant contribution, (d) The radiative corrections are reasonably impor- 
tant, comparable to the values oi_p\ coming from nonradiative multiple scattering. This is 
similar to the conclusion of Ref . [5| . 

An outline of the paper is as follows: In Sec. [Ill we briefly describe p^-broadening of 
a quark due to multiple inelastic scatterings and remind the reader that these inelastic 
scatterings can be recast in terms of the elastic scattering of a dipole in the medium. In 
Sec. lIIII we give, but do not derive, the general formalism which we use, a generalization of the 
formalism that Zakharov developed for radiative energy loss. In Sec. IIVI we quantitatively 
delineate the different kinematic regions for the radiated gluon. In Sec. |V] we evaluate the 
double logarithmic contribution to the average transverse momentum broadening, (p^)., of 
the quark due to gluon radiation. In Sec. |Vl]we evaluate the single logarithmic contributions 
to {p\). In Sec. IVIII we give the final result and do a rough estimate of the radiative 
correction. In Sec. IVIIII we discuss resummation and running coupling effects. In Sec. IIXI 
we give a derivation of our formalism and the relationship to that of Ref. (5| . 

II. p^-BROADENING AND DIPOLE SCATTERING 

A. A high-energy quark entering a nucleus from outside 

We begin by considering a high-energy quark passing through a large nucleus and scat- 
tering either inelastically, one gluon exchange, or elastically, two gluon exchange, with the 
nucleons in the nucleus. For the moment we work in the McLerran-Venugopalan (MV) |l4| 
approximation where there is no small- a; evolution. If the quark enters the nucleus with no 
transverse momentum and exits having transverse momentum px, then the distribution of 
transverse momenta is given by 

' ^e-'P^-^^Six^) (1) 



d^P± J (27r 

where S{x±) is the S'-matrix for elastic scattering of a quark- ant iquark dipole of transverse 
size x± on the nucleus. In the MV model the nucleons in the nucleus are uncorrelated so 
that 

S{x±) = exp [S'i(x_l)] (2) 

where Si{x±) is the elastic dipole-nucleon scattering S'-matrix which can be evaluated to be 

Siix^) = -^ (3) 

with the quark saturation momentum given by [2| 

Ql = ^^I^^LpxG{x,l/xl). (4) 



c 



In Eq. (jlj) L is the length of nuclear matter that the quark passes through, p the nucleon 
density, and G{x, l/x'jj the gluon distribution of the nucleon. While Eq. ([1]) is not a general 
formula it was shown in Refs. [3|, |4] that it is correct in leading and next to leading orders 
of small- a; evolution, thus carrying its validity well beyond the MV model. 



B. A high-energy quark produced in a nucleus 

The problem of p_L-broadening of a quark which is produced in a nucleus, say in a proton- 
nucleus or a nucleus-nucleus collisions, is quite different from that of a quark entering the 
nucleus from outside. For a high-energy quark entering the nucleus from outside, small-a; 
evolution corresponds to an evolved wave function of the corresponding dipole as it enters the 
nucleus. That is, small-x evolution mainly occurs before the quark, or corresponding dipole, 
reaches the target. However, a quark produced in the nucleus through a hard collision is 
produced as a bare quantum and small-x evolution is much more limited, and that evolution 
occurs in the nuclear environment. It is precisely these (small-x) radiative corrections to 
p_L-broadening of a (bare) quark produced in a hard collision which is the topic of this paper. 
At the MV level there is no distinction between a quark entering a nucleus from outside 
and a quark created in the nucleus so that Eqs. ([I]), ([2]), (E]) and (jl]) are valid while Eq. ([T]) 
should remain valid even in the presence of radiation in the nuclear medium, that is beyond 
the MV approximation. 



III. p^-BROADENING DUE TO GLUON EMISSION; THE GENERAL FORMAL- 
ISM 

We have seen that quark p^-broadening due to multiple scattering of the quark in a 
nuclear medium can be succinctly expressed in terms of the elastic scattering of a quark- 
antiquark dipole in the medium. The quark part of the dipole corresponding to the quark 
in the original amplitude and the antiquark in the complex conjugate amplitude. This same 
picture holds true when radiative corrections are included |3|, |j]. p^-broadening due to 
multiple scattering and single gluon emission is given by elastic scattering of a dipole in the 
medium where the dipole emits, and later absorbs, a gluon as it passes through the medium. 
We shall focus on soft gluon emission in which case the quark and antiquark making up 
the dipole can be viewed as fixed in transverse coordinate space during the emission and 
absorption of the gluon. The gluon, however, moves in transverse coordinate space during 
its lifetime, starting from either the quark or antiquark and ending again on one of the two 
parts of the quark-antiquark dipole. During the time that the gluon fluctuation lives, the 
quark- ant iquark-gluon system elastically scatters off nucleons in the nucleus. While this 
three body scattering can be done exactly in this section we shall limit our discussion to 
the large Nc and soft gluon limits where the picture is especially simple. A more general 
discussion is given in Appendix \M 
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FIG. 1. A quark-antiquark-gluon system moves in a nuclear medium with length L. The transverse 
coordinates of the quark and antiquark are and x± respectively. A gluon propagates from 
(i?i_i_,zi) to {B±,z), and its motion is described by the Green function G{B±, z; Bi±, zi). 



The essential part of the calculation involves the propagation of the quark-antiquark- 
gluon system through the nuclear medium. We choose the position of the quark to be at the 
origin in transverse coordinate space and the antiquark to be at x± and these positions are 
fixed during the evolution. Let G{B^, z; Bi±, zi) be the Green function for the propagation 
of the system from longitudinal position Zi to z where the gluon has transverse position 
Bij_ at Zi and transverse position Bj_ at z. (Since our problem is of relativistic motion 
longitudinal coordinates and times can be identified.) We normalize G by 



G{B^,Zi;B^^,z^) = 6^'\B^-B 



1-L> 



(5) 



G is pictured in Fig. [TJ In the large N^ limit the evolution of G is straightforward and is 
^iven by Ref . |5| , with u the energy of the gluon, 
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where q is the quark transport coefficient, q can be obtained from Eq. i^ a.s q = Q'j./L 
and the relationship to the gluon transport coefficient is ggiuon = qNc/Cp- The first term 
on the right hand side of Eq. ([6]) represents free gluon propagation while the second and 
third terms correspond to scattering of the quark- (antiquark part of the gluon) dipole and 
the (quark part of the gluon)-antiquark dipole in the medium. It is convenient to rewrite 
Eq. ([SD as 
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Eq. ([7]) is the equation for a two-dimensional harmonic oscillator with an imaginary 
potential. The solution, with initial condition Eq. (E]), is 
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In the large Nc limit the ojq in Eq. (fTOj) is exactly as in Ref. [5|. One gets a formula for 
the radiative contribution to 5'(x^), in Eq. ([1]), by connecting the gluon at ^i, and z to the 
quark and antiquark and by using the usual MV formula during the times when the gluon 
is not present. Suppose the nuclear matter is located in the longitudinal region Q < z < L. 
Then with 

S{x^)= f^N{x^,co), (11) 



N{x±,co) is given by 
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where Gq is obtained from G by setting wq = ^ = in Eq. (jHj). A few comments on Eq. (fT2l) : 
(i) The overall normalization will be fixed in the next section when we evaluate the single 
scattering term in G and connect the result with that from the standard dipole scattering 
formalism, (ii) The contribution in Eq. (1121) will give the radiative contribution to p±- 
broadening. One should add this contribution to the purely multiple scattering contribution 
discussed in Sec. [TTl (iii) The subtraction of the Go term in Eq. flT2|) subtracts the medium 
independent contribution, (iv) Eq. flT2l) is a generalization of an analogous formula intro- 
duced by Zakharov [8|, |9j to evaluate quark energy loss due to gluon radiation in a nuclear 
or thermal QCD medium, (v) Finally, a more complete derivation of Eq. flT2l) is given in 
Sec. [IX] and Appendix \^ 



IV. THE VARIOUS KINEMATIC REGIONS 

Before going on to evaluate the single and double logarithmic contributions to Eq. (TTTj) 
it will be helpful to qualitatively determine the kinematic regions in u and t = Z2 — zi which 
are important in Eq. (ITT!) and Eq. (fT2|) . Clearly < t < L, however when t < Iq we shall 
see that special care must be taken so that a straightforward use of Eq. fllip only occurs 
in Zq < ^ < -^- At Zi and Z2, Bij_ and B2± are taken to be either or x±. However, at 
2;- values between zi and Z2 one expects that 5_l may be much larger than Xj_, which we 
always suppose to be of size x± ~ ^/Qs- We may estimate the typical size of B± using free 
gluon diffusion between zi and Z2. This gives 

Bl ~ -. (13) 

However, when Bj_ gets too large, G, given by Eq. (|8]), will go to zero exponentially in 
B\. It is not hard to see that G will be small unless 

M<1 or t<y|. (14) 

Finally when Bj_ becomes as small as x\ the gluon will become very coherent with a part 
or all of the quark-antiquark dipole and, again, special care will be required. In terms of t 
and u the region B'j_ ~ x\_ corresponds to, using Eq. (fT3!) above and x^ ~ l/Ql ~ ^/(jL, 

The various boundaries and regions are shown in Fig. [2J The line (a) in that figure is 
given by t = cu/qL as in Eq. flTSj) . Line (6) is given by t = \/co/q as in Eq. flT^ while line 
(c) given hj t = Iq is given by the size of the nucleons in cold matter or by the Debye length 
in hot matter. 




FIG. 2. The phase space of the radiated gluon. t and co are the hfetime and energy of the gluon, 
respectively. The double logarithmic region is enclosed by the lines (a), (6) and (c). Regions Ai 
and A2 combined together are the region A. 

As we shall shortly see the region A, enclosed by the boundary lines (a), (6) and (c), is 
a region where Eq. flT2|) gives a double logarithmic contribution to Eq. flTTl) . This double 
logarithmic region is dominated by the single scattering contribution to N{x±,u). The 
region B, above line (6), is a region of strong absorption of the quark-antiquark-gluon system 
due to the separation of the gluon from the quark-antiquark pair being very large. In the 
next section we shall first evaluate Eq. ( TTTj) in the region A, a relatively easy task, to 
obtain the double logarithmic contribution. Then we shall, in turn, integrate Eq. flT^ 
across the boundary lines (a), (6) and (c) which, after doing the remaining logarithmic 
integration, will give the complete single logarithmic contributions. The constant, non- 
logarithmic contributions coming near the intersections of the lines [(a), (6)], [(a), (c)] and 
[(6), (c)] are beyond what we are currently able to do. 



V. THE DOUBLE LOGARITHMIC CONTRIBUTION TO (p^) 

In this section we shall evaluate the leading radiative contribution to 

dN 



{pI) = 

Using Eq. ([1]) it is straightforward to get 
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{pD = -^l,S{x^) 



(17) 
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which we shall now evaluate using Eq. f lTT]) and Eq. flT^ . In the double logarithmic region, 
region A of Fig. [H the typical values of B_i are much larger than Xj_ so that we evaluate the 
differences in Bij_ and B2±, taken at and x±, by expanding G so that Eq. flT^ becomes 
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(18) 



We have dropped the x^ term in the exponent since the x^ factor in Eq. (fT8|) is what is 
necessary to evaluate Eq. (fT7j) . It is straightforward to evaluate the derivatives of Bi± and 
Sax in Eq. ([18]) as 



(Vb,^ ■ Vb,S{g{B2±, t- Bi^, 0) - Go(S2x, t; Six, 0)) 




-Bix— -B2X— 

[4 - sin^ uot] - 4 \ (19) 



where the final term in Eq. flTUl) comes from Gq- Now taking u^t small one gets 



Using Eq. (EO]) in Eq. ([18]) along with Eq. ^ and Eq. ([I?]) 



~ TTt ■ 



(20) 



where the limits of integration in Eq. ( 1211) are given by the boundaries of the region A in 
Fig. [21 The result in Eq. ( |2T1) does not assume the large iVc limit. 

Finally, let us see that the answer Eq. ( 12T]) can be obtained almost effortlessly using 
the standard dipole picture for S{x±). We first note that it is natural to take x± < l/Qs 
otherwise S{x±) will be very small just from the elastic quark-antiquark dipole scattering in 
the medium. When x± ~ 1/Qs there is typically one quark-antiquark dipole scattering with 
the whole medium. Now in order to get a double logarithmic term from radiation one of the 
logarithmic integration must correspond to the various lifetimes of the gluon fluctuation. 
This means that the gluon lifetime is much less than the length, L, of the medium. Thus 
for there to be a significant chance for a scattering of the quark-antiquark-gluon system the 
gluon separation from the quark and antiquark must be very large to enhance the scattering. 
In this circumstance the single scattering, but only the single scattering contribution, can 
give a dBj_/Bj_ ~ dt/t logarithmic contribution. (Higher multiple scatterings would give 
extra powers of Bj_ thus destroy the logarithmic integration.) 
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FIG. 3. Single scattering in the dipole language. The emitted gluon with transverse coordinate 
B± is denoted by a double line in the large Nc limit. 

The single scattering contributions are shown in Fig. [3] for which the formula reads 



The first part of Eq. ( 122|) is tlie standard dipole probability distribution wliile tlie 

factor gives the (absorptive) single scattering amplitude for the quark-antiquark-gluon sys- 
tem to interact with some nucleon at impact parameter B±. In the double logarithmic region 
where B\ ^ x]_, Eq. fl22|) becomes 



X 



«,iv_. r^'^'^dB^ r^'^^^'^duj 



which, using Eq. (ITTl) . leads to Eq. (12T]) . The limits on the integrals come from the following: 
uj > lo/Bj_ is just the statement that the lifetime of the fluctuation, ujBj_, is greater than 
the nucleon size Iq. u < l/qB^ is the same as 1/B\_ > qB'^u, which is the transerverse 
momentum squared picked up by the gluon over its lifetime. Here we are requiring B± to 
be small enough so that multiple scattering is not important. The lower limit on the Bj^- 
integration is requiring that i?^ > x^ while the upper limit is necessary in order to have a 
nonzero range for the w-integral. The equivalence of Eq. ( l22l) with Eq. ( ITTl) and Eq. ( IT2l) at 
the double logarithmic level confirms the constant factor we have taken in front of the right 
hand side of Eq. (^^. 



VI. THE SINGLE LOGARITHMIC CONTRIBUTION FROM THE DIFFERENT 
BOUNDARIES 

In order to obtain a complete expression for (p^), apart from the double logarithmic 
term, we also have to calculate the single logarithmic correction to {p]_), which requires to 
do the dt and du integration in Eq. (12T1) over the entire allowed phase space, i.e. < t < L 
and qIq < u < qL"^. It is natural to divide the energy domain into two different regions, 
qll < u < qloL and qloL < u < qL"^, and do the calculation for each region separately. The 
reason we make this division is that the calculation near the boundaries (a) and (c) are very 
different. Moreover, when we do the du and dt integration over a certain region of the phase 
space, it is convenient to do the t-integral first but leave the u variable unintegrated. The tu- 
integration will be carried out after we combine contributions from different regions of phase 
space. The final expression for {p^}, Eq. ( l42l) . contains a double logarithmic term, as found 
in Eq. ( 12T|) . and a single logarithmic term, which comes from crossing different boundaries 
(a), (6) and (c), and an undetermined constant term. The following calculation consists of 
three different parts. First, we cross the boundary (6) by starting from some arbitrary point 
to in region A, then do the t-integral from to to L, which allows us to collect the double 
logarithmic contribution from the region A as well as the single logarithmic contribution 
from the boundary (b). Next, we deal with the boundaries (a) and (c) separately. The 
t-integration for those two boundaries will be done from to to such that when we add the 
contributions from the three regions together there is no to dependence. Finally, we do the 
w-integration and obtain our final expression for {pY). Note that crossing the boundary (6) 
is relatively simple and in our approximation Eq. ( 1T8l) is the correct formula that we should 
use in the regions A and B. So in the following calculation we first calculate the contribution 
to {p\) from regions A and B, then extract the single logarithmic contributions from the 
boundaries (a) and (c) by modifying the dipole formula Eq. (|22|) . One reason we have to 

10 



reformulate the dipole formula is that in Eq. (1T8|) a time integral is used, so when crossing the 
boundaries (a) and (c) we have to express the dipole formula in terms of a time integration. 
After all if we want to connect calculations from different regions of phase space, we have to 
connect them by the same variable. The complication arises from converting the transverse 
coordinates integration in the dipole formula to a time integration. Another reason is that 
below those two boundaries certain kinematic limitation of the dipole formula will be met, 
the formulas no longer take the simple and compact form. We will discuss the limitations 
and modifications in more detail as we carry out the calculation. 



A. Crossing the boundary (6) 



As we have seen in our previous calculation, Eq. flT8|l in the limit uot <^ 1 gives the same 
double logarithmic contribution as the single scattering contribution from the simple dipole 
calculation, Eq. fl22l) . The line t = \/co/q in Fig. [2] can be taken as the separation of the 
single scattering and the multiple scatterings. As we start from the region A, crossing the 
boundary (6), we enter the region B where the multiple scatterings play the major role. 
In the multiple scattering region the Green function is exponentially suppressed in Bj_ as t 
approaches L, so there is no double logarithmic contribution from this region. We can choose 
some initial point to, which is in the double logarithmic region below the boundary (6) but 
above the boundary (a) or (c) and satisfies wq^o ^ 1- Since now we are not only interested 
in the double logarithmic contribution from the region A but also the single logarithmic 
contribution from the boundary (b), we should use the complete expression Eq. (TT9l) in 
Eq. (TT8|) to do the t-integration from to to L. That is 
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It is straightforward to do the time integration, and Eq. fl24j) becomes 
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With Eq. fl25|) and using Eq. flTTl) and Eq. flT7|) . we can calculate {p\), but leave the cu 
variable unintegrated, that is 



u 



d 

duj 



(p1> 



-VlN(x^,oj) 



boundary b 



aMr. 



vr 



qL 



In 



(26) 



which not only contains the information about the double logarithmic contribution from the 
region A, but also the single logarithmic contribution from the boundary (6). 

B. Crossing the boundary (a) 



Our starting point is the dipole formula, Eq. fl22|) . but with some modifications. Near 
the boundary (a) the transverse size of the gluon is comparable to that of the dipole. Note 
that in Eq. ( 122|) a singularity occurs if we take |-Bx| -^ \^i-\- This is because we did not 
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include the usual virtual dipole scattering term, guaranteeing probability conservation, in 
Eq. ( 122|) . Indeed, the virtual term does not give a double logarithmic contribution and could 
be neglected in arriving at Eq. ( 123|) . However, the virtual term is important here and can 
be put in simply by replacing the —\qL\B\ + {B^ — x^)^] term in Eq. (1221) by 



f[5i + (5, 



XA_f 



X 



}-^B. 



[B± - xx) 



(27) 



a result which also follows directly from Eq. flT2|) . Moreover, we have to convert Eq. fl22|) to 
an expression with an integral over time so that the same to appears in order to match onto 
Eq. ( 12B|) . Note that x']_/B'j_{B± — x±Y is the gluon emission probability and is calculated 
from squaring the wave-function that a dipole emits a gluon. So we can keep the time 
integration in the wave function and write the gluon emission probability as 



Xi 



Bl{B^-xJ^ 



"T 



dti f^ dU 

tl JO ^2 



'^-^Blu/2U^^ _ g-i(i?^-.^)V2ti(5^ _ ^^) 



'^^u.Bl|2t,^^ _ ^iu(B^-.^f/2t,^^^ _ ^^^ 



(28) 



Substitute the modifications, Eq. (|27|) and Eq. fl28|) . into Eq. f l22|) . do the d'^B^ integration 
first, then introduce t = t2 — ti and do the ti integration from to —t. The result is 



uj-—S{x±] 
du) 



asNc 
' 127r 



gLRe 



to 



dt 



x\ iu}{x\Y 
T^ 2t2 



^iujx^/2t 



Ai 



+ -(1- 



,jiu}x^./2t 



UJ 



(29) 



The integral in Eq. f l29p goes up to to iii order to match onto Eq. (1261) . It is straightforward 
to do the t-integration and then take the limit where to is relatively large, one finds 



d asNcx\ 

uj—b[x^) = qL 

CO Atx 



In Mo + 7 - o 



(30) 



where uq = —iux\/2tQ and 7 is the Euler constant. The detail of the calculation is shown 
in Appendix IB 1[ Substituting Eq. f l30|) in Eq. flTTl) . we obtain 



-^<^^^> 



a.N, 



boundary a 



TT 



-qL 



In — T + 77 



bJX , 



7 



(31) 



Following the prescription in Refs. [21, ll3| and the discussion in Sec. Uwe have taken the 
V^^ in Eq. ( 1T71) to apply only to the x\ prefactor on the right hand side of Eq. (l30l) and 
not to the Inuo term. With this prescription, and after taking an appropriate scale for x\ 
instead of setting x\ = Q, {p\) in Eq. fl3Tl) corresponds to a typical p]_ and not the mean p\ 
which is not a useful concept in the present context. 

C. Crossing boundary (c) 

In order to see the limitation of Eq. fl22|) near the boundary (c) let us first determine 
the x-value of the nucleon gluon distribution probed by the high energy quark with a soft 
gluon radiation. Suppose the nucleon is at rest, P^ = (m, 0, 0, 0) where m is the rest 
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(a) (b) (c) 

FIG. 4. Induced gluon radiation of a high energy quark by a single scattering. 



mass of the nucleoli, and the quark is moving along the 2;-direction with a four-momentum 
Pm ~ (P)0,0,p). Consider a typical scattering diagram, for example. Fig. 4(a) /_ from the 



nucleon is very close to the value k\/2kj^. Then the x-value of the gluon distribution of the 
nucleon is 

The lifetime of the gluon fluctuation is t ~ 2uj/k']_, and we can roughly take m ^ l//o- 
Using the above approximations, we see that x ~ lo/t. That is near the boundary (c) the 
X- value of the gluon distribution of the nucleon approaches 1. It means that the transport 
coefficient q, depending on the gluon distribution, can no longer be taken as a constant and 
its x-dependence cannot be neglected in the calculation. 

Here we evaluate (p^) of the quark by identifying it with the average k'j_ of the radiated 
gluon. It is sufficient to take only the real diagrams into account, as shown in Fig. HI One can 
show that the virtual diagrams do not contribute to {p\) in the approximation |/_l| > \k±\. 
It is convenient to do the calculation in momentum space. Therefore from Fig. H] we have 



dN 



aMr. 



u 



dijj(Pk I 



vr^ 



pL 



dHi 



I' 



da 



ki{k^-i±_ydnj 



(33) 



where a is the high-energy quark-nucleon scattering cross section coming from the elastic 
scattering and the factor {asNc/7i^)l^/kj_{k± — l±)^ is the probability that a high energy 
quark emits a soft gluon. From Eq. (fT6!l we have 



(pi> 



boundary c 






a.Nr 



TT^ 



'-pL 



da 



—dH 
u ^dH 



d^k\ 



kj_lj_ 



klik^ - l±) 



2' 



(34) 



Let us write the scattering probability in terms of time integrals in order to separate the k± 
and l± dependence. Later the ^^-integration will be converted into a x-integration, and the 
/^-integral will give the gluon distribution of the nucleon. So we write 



2 72 



d^ki 



kU 



1 ^0 /"oo 

i- / dtj dh kld'k^ \k^e''"^'^/^^ - (k^ - /^)e^(fe^-'x)^*i/2-l 

^ J-oo Jo 

(35) 



fc2_(fcx-/±)2 4^2 

• [k^e-''"^''/'^ - {k^ - /^)e-^('=x-'x)'*2/2.] _ 
Introduce t = t2 — ti and do the ti-integral from — t to 0, Eq. ( 135|) becomes 

d^k 



/ 



k' 



.h _ g-ifcito/2c^^ 



(36) 
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where we have dropped e'^-L*''/^'^ factors not enhanced by sniall-A;^, integration. The de- 
tail derivation of Eq. ( I36l) is given in the Appendix IB 21 Now we see that in Eq. (l36ll 
a fc^-integration appears. It can be converted into an integration over x using Eq. (1321) . 
Furthermore, the /^ factor can be used to obtain the gluon distribution of the nucleon in 
Eq. ([34]) via 



£hll 



da in ttsCp 



xQ 



(37) 



where xQ contains the general x-dependence of the gluon distribution of the nucleon. Putting 
the above two steps together we obtain 



(pi> 



a.N, 



boundary c 



TT 



^ f dcj f dx , 

pL — — (1- 

J (^ Jo X 



-irnto 



m-1 



xQ. 



(38) 



We shall separate the gluon distribution into two different parts, one part is small at 
small-x while the other one is the gluon distribution used in the transport coefficient and 
has the usual small-x dependence. So we write 

xQ = [xQ - xG) + xG (39) 

where xG is the usual, small-x, form of the gluon distribution as used in Eq. (j4]). Using 

1 



Re 



I dx^^ 
Jo ^ 



-ixmto 



) = 7 + ln(mt, 



oj 



O 



mt, 



(40) 



with to being large, we find 
d 



-i^"^) 



''-gL[ln(mto) +7] + 



boundary c 



TT 



TT 



^-r i^dx 

qL / — 







X 



xQ 



(41) 



VII. FINAL EXPRESSION FOR 



(pl) 



Now we shall gather all our results from Sec. IVI A[ Sec. IVl Bl and Sec. IVICI In the region 
qIqL < uj < qL"^ the t-integration encounters the boundaries (a) and (6), so we add Eq. (EI]) 
and Eq. ( 126]) and do the w-integration from qloL to qL"^. In the region ql^ < uj < qloL the 
t-integration encounters the boundaries (c) and (6), so we add Eq. ( 14T]) and Eq. (p6]) and 
do the w-integration from ql^ to qloL. Finally, we can combine these two results together 
and obtain our final expression for the radiative corrections to the transverse momentum 
broadening 



(A) 



— — qL In -^ + 



TT 



qL 



In 



SttlIq 
xP'qL 



1 

3 + 



^ dx 

X 



xG 



In^ + C 



(42) 



where the constant G cannot be determined in our current calculation. We remind the 
reader that qLx^/A in the second term should always be taken to be of order one. (See the 
discussion in Sec. HI) Note that if we rescale lo, i.e. lo ^ c ■ lo, where c is some arbitrary 
constant, then the first two terms in Eq. (jl2]) remain the same, however, the constant G is 
changed by an arbitrary amount. Therefore we can only determine {p\) up to some arbitrary 
constant, but the double and single logarithmic terms in Eq. (j^2l) are exact. 
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In order to do a rough estimate of the size of the radiative correction to p^-broadening of 
a quark in hot matter we take L ~ 5 fm, m ~ 300 MeV and m/o ~ 1 as well as xQ jxG ~ 1. 
(For cold matter xQ jxG would be much smaller, perhaps something like (1 — x) 
with x\qL = x']_Q'l ~ 4 and a^ ~ 1/3 



Then 



(A) 



rad 



0.75 qL, 



(43) 



a rather large contribution from radiative corrections. The dominant part of Eq. (1431) comes 
from the double logarithmic part and is clearly sensitive to the exact values one takes for 
as, L and Iq. 



VIII. RESUMMATION AND RUNNING COUPLING 

In this section we shall briefly discuss two extensions (modifications) of our analysis. 
The first of these is a resummation of our calculation while the second touches a somewhat 
different calculation because of the way in which we shall approach running coupling effects. 

A. Resummation 




ql„ qL ^L 

FIG. 5. The double logarithmic region and its boundaries in (x,A;^) rather than (t,tj) of Fig. [2j 



It is not difficult to resum all leading double logarithms. To do this it is convenient to 
draw Fig. [2] in a somewhat different way using variables familiar from double logarithmic 
resummation in deep inelastic scattering, namely k^ and x as defined in Sec. IVICI The 
boundaries of Fig. [2] now appear as in Fig. [51 The boundary (6) now is given by 



X 



qk 



Ql 



(44) 



where, in the double logarithmic approximation, we do not distinguish the mass m in Eq. (I32ll 
from Iq. We shall begin by repeating the double logarithmic calculation given in Sec. IVlin 
terms of our current momentum variables. We may take the expression for {p^} directly 
from Eq. flMj) which, using Eq. (^ and Eq. ( 137]) with xQ replaced by xG, can be written as 



/ 2\ asNc^^ Z"^' dk"^, 
[P±) = -^qL 



TT 



Qt 



B 



dx 

l/kl X 



(45) 



15 



and leads to 



(pi> 



471 ^ 2! V ^0 



(46) 



exactly as in Eq. fl2T]) . The integration region in Eq. fH5|) goes over the region A in Fig. |5l 
Graphically it is convenient to view Eq. fH5|) as corresponding to the graph shown in Fig. 6(a) 



where l'j_/k'j_ is taken as usual in the logarithmic approximation. From the point of view 
of how the Z-lines hook onto the quark- ant iquark pair the quark-antiquark pair acts like a 
gluon because the quark-antiquark separation, x±, is much less than l/l±. The next term 





(a) (b) 

FIG. 6. Schematic illustrations of Eq. ([^5]) and Eq. ([T7|) respectively. 



in the leading double logarithmic series is illustrated in Fig. 6(b) and is given by 

•1 _7__ PX^ J.. / _ AT \ 2 



(p1> = qL 



asNc 



vr 



<3» dkl^ 



'''^^ dklj_ 



k'i 



dxi 



i± JQl/kl^ ^1 



'(ll^2A 



dX2 
X2 



qL 



OsNc 

Air 



1 , 4 (L' 



The pattern is now clear and the full leading double logarithmic sum is given by 



{pD = qL 



Arc 



asN, 



-h 



In 7^ 

'■n 



asN, 



TT 






(47) 



(48) 



where the nonradiative, qL, term is included in Eq. fHS]) and Ii{x) is the modified Bessel func- 
tion. The effect of the boundary (6) is seen by noting that if the ordered x and k'j_ integrals 
went between Iq/L < x < 1 and Ql < k'j_ < Ql without the restriction on the lower limits of 
the X integrals in Eq. fllS]) and Eq. (1171) the answer would have been qL Iq[-\/ asN^ Tiln^L"^ /11)] 
15| rather than the right hand side of Eq. (jlSj) . Though the resummation is interesting in 
its own right, for reasonable values of the parameters, say L/Iq ~ 5, a^ = 1/3 the correction 
Eq. (H6|) is dominant and the resummation does not appear useful. 



B. Running coupling 



In all that we have done so far we have taken q to be scale independent and we have 
not included any scale in the coupling a^. Naturally, the scale in the as in Eq. ( H2l) should 
be Q^. However the scale in q in Eq. fH2|) is not straightforward. For example to solve the 
region near boundary (6) (see Fig. |2]) we would have to solve Eq. (jTj) for a scale dependent 
q to get a generalization of Eq. ([8]) and that appears beyond what we are able to do. We 
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can, however, evaluate the double logarithmic contribution using a running coupling, and 
the result has same interesting features. To include running coupling effects at the double 
logarithmic level we take the q in Eq. (H5|) and evaluate it at a scale k^ inside the integral 
over k\. Even this however is not completely unambiguous. In Eq. (jl]), using Ql = qL we 
see that the scale comes in both in a^ and in xG. We will keep the scale dependence in a^ 
but drop it in xG. (In this context we note that the scale in xG would involve a lnln(Qg/A^) 
while as{Q^) would be given by as{Q^) = l/(61n(5s/A^) so that the scale dependence in xG 
is rather weak.) Then Eq. fH5|) becomes 



/ 2\ «.^c., Q^ f^^dkl 1 f' dx .^^. 

(P±) = -^^I^tI / .. TT-l^ /„.„-' (49) 



^ ^' -iQl kl In K Jq2/,2 X 



A2 



where now the q in Eq. ( H9|) is evaluated at scale Q^. Identifying A^ and Ql at the double 
logarithmic level one finds 

What is rather remarkable here is that there are no In In terms coming from the fc^- 
integration. The leading double logarithmic result here is simply a factor 2 larger than 
in the fixed coupling case. Probably we should take this as a measure of uncertainty for our 
whole calculation. Higher orders in {asNc/Tc)\n^{L'^/ll) can be evaluated but do not seem 
to lead to any simple result as happened in Eq. (HHj) for fixed coupling. 



IX. DERIVATION OF OUR FORMALISM 



{xi±,Zi) {X2±,Z2) 



J_ 



< \ x± 



I I 

' /x/X/^ {B2^,Z2) 



(Six, zi) r\y^ I 
\ > \ 0^ 

oooooooo 

L 

FIG. 7. A quark-antiquark-gluon system in the case when the emitted gluon carries a finite fraction 
of the quark energy. 

In this section we outline the main steps of deriving Eq. (TT2|) in the case where the emitted 
gluon, real or virtual, carries a finite fraction ^ = u/E of the quark energy. More details 
of the derivation can be found in Appendix |Al Let us start with the situation where the 
gluon is already present in the medium, as illustrated in Fig. [71 At z < Zi and z > Z2 
the transverse coordinates of the quark and antiquark are fixed. The contribution of the 
multiple scatterings of the quark-antiquark pair can be written as the product of two dipoles 
of transverse sizes xi_l = x±{zi) and X2± = Xj_{z2) multiply scattered by the medium. This 
dipole type multiple scattering contribution becomes the phase factor in the final expression 
Eq. fl55|) . During the lifetime of the gluon, i.e. Zi < z < Z2, the quark-antiquark-gluon 
system is multiply scattered by the medium, and the transverse coordinate of the center of 

17 



mass of the quark-gluon system is fixed relative to that of the antiquark. As a result one 
has 

X2± - ^B21_ = Xi_L - ^Bij_. (51) 

The propagation of the quark-gluon system in the medium can be described by a Green 
function, G{B2_l, Z2\ Bi±, zi), which can be found from the Hamiltoniaio 



^ ^ -mtr) - iH^- ^ ^[^- ^ '^^^ - ^=^' 



2 2 1 



(52) 



In fact, G{B2±, Z2] -Bi±, Zi) does not depend on X2±; the X2_L-dependence should be eliminated 
according to Eq. ( 15T]) . Moreover, the emission and absorption of the gluon from the quark 
or the antiquark can be taken into account by applying the operator 



^g^Cp 



^.(0(Vi.,x - ev.,J ■ (Vb,, - evx,j 



e(i-0' 

to the Green function with the splitting function given by 



(53) 



n(o 



e 



(54) 



Putting the three contributions together, we find 



S{x2^) = ^Reldi ^^(^) 



E' J ^ai-0 .0 



L i>Z2 

dz2 I dzi 



X 



g-i,>?,.,+.i,(L-..)]y^^^ ■ (Vb,, -^V,,0G(52±,;22;5i±,^i) 



Vb2_l ■ (Vsi^ - i'S7 xi^)Gq{B2u Z2] Bi^, zi] 



■B2±=3:^2±,-Bi±=0,a;ix = (l-?)a;2± 



B2^=0,Bij^=0,xi^=X2 



(55) 



Note that -Bi±, -B2±, xi±_ and X2± obey Eq. flSTj) . The p_L-broadening of the quark is given 
by 



{pD = -yl,s{x2) 



(56) 



Xo I =0 



It is easy to see that in the soft gluon hmit, ^ — )■ 0, Eq. (|55|) reduces to Eq. (IT2|) . 
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^ In the soft gluon and large N^ limit the Green function defined in this section is actually the complex 
conjugate of that in Eq. ([6|). However, they give the same answer since only the real part of G is needed 
to calculate {p\)- 
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FIG. 8. Medium-induced gluon emission (Fig. 8(a)) and quark self-energy (Fig. |8(b)[ ). Both 
diagrams are composed of the multiple scattering of the quark-antiquark dipole at z > Z2 and 
z < zi and of the quark-antiquark-gluon system sX Z2 > z > zi. During Z2 > z > zi^ a, h and c 
respectively stand for the antiquark of energy E, the quark of energy (1 — ^)E and the gluon of 
energy w = £,E. 



Appendix A: Review on the path integral approach of Refs. 



In Zakharov's formalism [8|,|9| (see also Ref. 16| for a discussion similar to the one below), 
the propagator of a parton 2 in a nuclear medium takes the form 



Gi{x2L,Z2]Xi^,Zi) 



X ^{z\)=Xl^, X ^{Z2) = X21_ 



Vx± exp <^i dz 



E; 



«^2 



x, - VAx^^z) 



(Al) 



where Ei is the energy of the parton and Vi{x±,z) is the potential of the medium. The 
momentum eigenstate of the parton in the vacuum is 



zi 



il)ip^{x±,z) = i{x±,z\pj_) = e^^'""^ *2E,^ 



(A2) 



We first derive explicitly the formula of S{x) for real gluon emission process (Fig. 8(a)) and 



for the quark self-energy (Fig. 8(b)) inside the medium. For both processes the emitted 



gluon carries a finite fraction ^ = co/E of the quark energy E. The Feynman rules for 
evaluating the diagrams in Fig. |8] are: 

• At 2; = and z = L, the wave functions of the quark and antiquark are given by 
Eq. (lA2p and its complex conjugate, respectively. 



The propagators of the quark and gluon are given by Eq. flAip and the propagator of 
the antiquark is given by the complex conjugate of that of the quark. 

The product of the two gluon emission/absorption vertices at z = Zi and z = Z2, after 
summing over the two polarization states of the gluon, is given by 

/ (fxi±(fx2±l{XbiJ_, Zi] Xcl±, Zi\Hl{Xi±, Zi) \Xal±, Zi) j l{Xa2±, Z2\Hl{x2±, Z2) \Xb2L, Z2] Xc2±, Z2) j 

cPpi± (Pki^ cPp2L (Pk2± 2g'^CF 



X e 



(27r)^ {2nr (2vr)^ (3.), ^(^ _ ^)^.(0(/^ix - ^m) " (A:^^ - ^m) 

ikl±-iXcl±-Xi,i±)~ipi±-(Xal±-Xbl±)-ik2±-{Xc2±-Xt2±)+iP2±-{Xa2±-Xb2±) 



(A3) 
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Integrate out all the transverse coordinates and the time zi and Z2. At the end the 

dN 
(PpA 



contribution from these two diagrams to -S^ is 



dN, 



± 



1 



''A{2T:mi-i)E^A 



d?p 

d^PL Jo 



1 



Ai27rna~OE^A 



X 



X 



Re(diagram in Fig. 8(a) with Z2 > Zi), (A4) 



Re(diagram in Fig. 8(b) with Z2 > zij, (A5) 



where A is the area of the transverse cross section of the medium. The distribution of the 
transverse momentum of the quark in the final state is given by 



dN _ dNrg dNyg 



d'^p± d'^p± d^pi 



(A6) 



Eqs. flA4p and flASp can be simplified significantly by taking the same medium average (■ ■ •) 
as that in the BDMPS formalism. In this case all the path integrals can be solved analytically 
in the saddle-point approximation. 

The S'-matrix for the multiple scattering of a quark-antiquark dipole is 

S'j (X2_L , X2_|_ , Z2:, Xi_L, X;^_|_, Zi) = {Gi{x2±, Z2', Xi±, Zi)G*{x2j_, ^2! 2;^_|_, Zi)) 



Vx±Vx'j_ exp <i dz 
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yxi -Vi{X:i,z) 



l>Z2 

I I dz' 

'zi 



^x'l-V,{x\,z') 



I 



Vx±Vx' ± exp < dz 



Z2 



Zl 



■ Ei .2 Ei ■ 2 1 I \ ( \ 

zyx^ - ^ Y A - 2^"(^-L ~ ^±' z)n[z) 



where 



CTi 



i(a;_L -x'^,z)n{z) = {\yi{x^,z) -Vi{x'^,z)\^J ~ -^ql^^ - x'^\ 



(A7) 



(A8) 



The boundary conditions x'^{zi) = x[j_, x'^{z2) = x'^^-, xx_{z\) = Xi± and x±{z2) = X2± are 
used in Eq. ( I A 71) . It is easy to show that 



(A9) 



f d''xi^d%^S,{x2±, x'2^, Z2; xi^, x[^, zi) = e-3«(-2-^i)l-2^--2xl' 
f d^X2±d^x'2^e-'P^<^^^-<^^S,{ 



X2±,x'2^,Z2;xi±,x[^,zi) = e-*P^-(^i^-<x)-49(^2--i)l^ii-4xl . 

(AlO) 



During Z2 > z > Zi one needs to evaluate the following path integral 

Gabc{Xa2±, Xb2±, Xc2±, ^2] a^al±, ^hX^^ a^clX, Z\) 
= {Gb{Xb2J_, Z2] Xbi_L, Zi)Gc{Xc2±, Z2] Xcl±, 2l)G'*(Xa2±, Z2] Xal±, ^l)) 



VxaiT^XhiVxci^ exp < dz 



-J^l^ + Y^c± - Y^^^ ~ H(a;fe±, z) - Vc{Xc±, z) + Va{Xa±, z) 



■ Ea .2 1 



i^^ll. + ^Y*^^ ~ *T^"^ ~ l^(^ahc[XahL,XcbL)n[z) 



(All) 
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where 



with 



^Crabc{Xab±,Xcb±)n{z) = {\VbiXb±,z) + Vc{Xc±,z) - VaiXa±, Z)\'^) 
- -^qI Xlb± + ^ [xlb± + {Xcb± - Xab±y - Xl^j_] \ (A12) 



Xab± = Xa± - Xb±, Xcb± = Xc± - Xb±. (Al3) 



Here indices a, b and c stand for the antiquark of energy E, the quark of energy (1 — C,)E 
and the gluon of energy u = C,E, respectively. The boundary conditions in the definition 

of Gabc in Eq. (|A11P are Xa±{Zi) = Xal±, Xa±{z2) = Xa2±, Xb±{Zi) = Xbl±, X b±{z2 ) = Xb2±, 

Xc±{zi) = Xci± and Xc_l{z2) = Xc2±- The evaluation of the path integral in Eq. (lAlip amounts 
to calculating the classical trajectory of the partons with the Lagrangian 

E . i 

L = —^(1 - Oxlb± + E{^Xcb± - iab±) ■ Xl_ + -(Tabc{Xab±, Xcb±)n{z), (Al4) 

where 



As a result one finds 10 



Gabc = Gl{^Xc2± + (1 - 0Xb2±, Z2] ^Xcll. + (1 - Cl^blL, Zi) 

^ Ga{Xa2Li Z2'-, Xal^-, Zi)G{Xcb2±-, ^2] a^cblX, ^l), (A16) 

where the propagator G is defined by the Hamiltonian 



2Eai-0 4 
with 

Xab±iz) -^Xcb±iz) 



yi Xlb± + ^ [xlb± + {Xcb± - Xab±)^ - Xl^j_] I (A17) 



Z — Zi 



and 



= Xab±{Zl) - C,Xcb±{Zl) H [Xab±{z2) " C,Xcb±{z2) - Xab±izi) + C.Xcb±{Zl)] 

Z2 — Zi 

= Xabl± — C^icbl. H \Xab2± — ^Xcb2± — Xabl± + C^cbl^ , (^18) 

Z2- Zi^ 



G°(xx, z; 0, 0) = ^e'^eiz). (A19) 
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Insert (lA9p . (lAlOp and (1A16P into (1A4I) and (lASp . and after some algebra we have 



X V 



dz2 / dzie-49[(i-5)^^i+(i'-^2)]xL 






a;cbl±=0, Xabl^={'^-i)x2^, Xcb2±=X2± 

(A20) 






^2 



X V.,,,, ■ (V 



l,n.-.,..2^ 



b2± V * a;ci,i± 



^'^XM±)G{Xcb2±, Z2] Xcbl±, Zi] 



a^c62J_— 0— XcblXi Xabl_L—Xab2A_—X21. 

(A21) 



and S'(x2±) = Srg{x2A.) + 'S'i,g(x2±). Here G is the propagator defined by H^c in Eq. flAlTp 
with 

XabL{z) = Xabl± +^{Xcb±{z) - Xcbl±) (A22) 

and as a result G only depends on Xabi±, Xcbi± and Xcb2±- Eq. (I55l) can be obtained by 
replacing Xa6± by x^ and Xcb± by i?_L and subtracting out the limit q ^ 0. 



Appendix B: Two integrals 
1. First one 

In this appendix we evaluate the integral 

h= [d^B^B^-{B^-x±) 



Xi 



Bl{B^-x^) 



(Bl) 



With the help of Eq. ( !28|) Ji can be written as 



h 



u 



2 /-o dti r^ dh 

t: . ^r^ ''1 .In 1^1 



-oo ''I Jo 

B± ■ {Bj_ — x±) exp 
1 



-r IT -i d'B^B^-{B^-x±){Biexp 



t-\ to 



Bleu , .{B±-x±yu 
"2t7 



+ 1 



+ iB± - x±) exp 



t-iBi-x^Yuj 



2t2 

1 1 



B± ■ {Bj_ — Xj_) exp 



Blu XB^-x^)^u' 



2U 



2ti 
(B2) 



The first term and the last term in Eq. flB2p give the same result after the integration. It is 
straightforward to do the -B^-integration, the first term in Eq. (]B2p gives 



4:Tri 



dti I dt2 ^\ 

h {t2 - hy 



Ani 



to 



dt 



UJ Jo t^ 



dtiti(ti +t) 



TCXj 
3Mn 



(B3) 
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where uq = —iujx]_/2tQ. The second and third term in Eq. (lB2p give the same result. Again, 
one can do the i^^-integration first and find 



TTO; 







dti i dt2 

-oo Jq 



oo ^iujx\/[2{t2-ti)] 



2 .i_^r\L ""^ {t2-tiY 



otit2 a;_|_ otit2 2 ■ ^1^2 4 

'u''{t2-h)~~~u{t2-hy''^~\t2-hy''^ 



■ (B4) 



Introduce the time difference t = t2 — ti and do the ti-integral from — t to as in Eq. (1B3I) . 
then Eq. (IB4p becomes 



n 



to 



^^gia;xi/2i 



— + ^ 



a; 2t 4t2 



(B5) 



With the help of the incomplete gamma function T{a, x), we can write the three t-integrals 



as 



to 



dt e^^^i/2i 



lUX , 



■r(-i,Mo), / dt 



to piuix\/2t 



fto J,uix\/2t 2i 

r(o,Mo), / dt — - — = — ^r(i,Mo)- 



t^ ujx\ 



Then Eq. (jBSP becomes 



TT 



r(-i,no) + ^r(o,Mo)-^r(i,'uo) 



Finally, combining Eq. ( JB6I) and Eq. ( JB3I) we arrive at 



/i 



2^x^ 



- - r(-i, no) + ir(o, no) - ^r(i, no) 

no 2 2 



(B6) 



(B7) 



Using the properties of the incomplete gamma function r(a + 1, x) = aT{a, x) + x"e~^ and 
the expansion r(0,a;) = —7 — Inx, for a; <^ 1, we can further simplify Eq. flBTD in the limit 
no — 7- as 

Ji = -TTx^ In no + 7 - - . (B8) 



2. Second one 



Here we evaluate the following integral 



Qi rC \ rv 



I' 



Again we rewrite it in terms of time integrals as in Eq. (135|) . and find 

/2 = __!_ r dt, r dt2 [d^k^kll A;^e^'=i(*i-*2)/2c. + (^^ _ ^^Y^^(k^-l^)Ht,-t,)/2u 

— k± ■ {k± — l±) exp 

— k± ■ {k± — l±) exp 



2u; V ^ ti-ta/ 


Ji tit2 
2ujti-t2_ 


2u; V ti-ts/ 


. l±_ tit2 
2u ti — ^2 



(B9) 
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For the third and forth terms in Eq. (lB9p . we shift the k± variable according to k± — t- 
kj_ — t2l±/(ti — ^2) and k± ^ k± + til±/(ti — ^2) respectively in order to make all the terms 
to have the same phase factor. Averaging the direction of k±, we find 



4c^ -/ -/-oo JO 



/"OO 



d'^k± I dti I dt'2 



\ 2{tl+t2f 



fcl + -^i^(A;l)2|>e-^'=i*/2-. (BIO) 

(jj ti — t2 



(tl - t2 

Introduce t = t2 — ti, and it is straightforward to do the ti-integration which leads to 



/2 = — 

60J 



where in the last step we have neglected the e'^^x*"/^"^ terms. 



[1 

[2 

[3] 
[4] 
[5 

[6; 

[7 

[8] 
[9] 

[lo; 

[11 

[12 

[13 
[14; 
[15 
[16 



G. T. Bodwin, S. J. Brodsky and G. P. Lepage, Phys. Rev. Lett. 47, 1799 (1981). 

R. Baier, Y. L. Dokshitzer, A. H. Mueller, S. Peigne and D. Schiff, Nucl. Phys. B 484, 265 

(1997) [hep-ph/9608322]. 

Y. V. Kovchegov and K. Tuchin, Phys. Rev. D 65, 074026 (2002) [hep-ph/01 11362]. 

A. H. Mueller and S. Munier, Nucl. Phys. A 893, 43 (2012) [arXiv: 1206. 1333 [hep-ph]]. 

B. Wu, JHEP 1110, 029 (2011) [arXiv:1102.0388 [hep-ph]]. 

R. Baier, Y. L. Dokshitzer, A. H. Mueller, S. Peigne and D. Schiff, Nucl. Phys. B 483, 291 

(1997) [hep-ph/9607355]; 

R. Baier, Y. L. Dokshitzer, A. H. Mueller and D. Schiff, Nucl. Phys. B 531, 403 (1998) 

[hep-ph/9804212]. 

B. G. Zakharov, JETP Lett. 63, 952 (1996) [hep-ph/9607440]. 

B. G. Zakharov, JETP Lett. 65, 615 (1997) [hep-ph/9704255]. 

R. Baier, D. Schiff and B. G. Zakharov, Ann. Rev. Nucl. Part. Sci. 50, 37 (2000) [hep- 

ph/0002198]. 

A. H. Mueller, Nucl. Phys. B 415, 373 (1994). 

Y. V. Kovchegov and E. .M. Levin, "Quantum Chromodynamics at High Energy" , Cambridge 

Univ. Press (2012). 

P. B. Arnold, Phys. Rev. D 80, 025004 (2009) [arXiv:0903.1081 [nucl-th]]. 

L. D. McLerran and R. Venugopalan, Phys. Rev. D 49, 2233 (1994) [hep-ph/9309289]. 



m- 



See Sec. 2.4.6 in 

J. -P. Blaizot, F. Dominguez, E. lancu and Y. Mehtar-Tani, JHEP 1301, 143 (2013) 

[arXiv:1209.4585 [hep-ph]]. 



24 



